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BORN-OPPENHEIMER APPROXIMATION FOR AN ATOM IN 
CONSTANT MAGNETIC FIELDS 

SOHEI ASHIDA 


Abstract. We obtain a reduction scheme for the study of the quantum evo¬ 
lution of an atom in constant magnetic fields using the method developed by 
Martinez, Nenciu and Sordoni based on the construction of almost invariant 
subspace. In Martinez-Sordoni m such a case is also studied but their re¬ 
duced Hamiltonian includes the vector potential terms. In this paper, using 
the center of mass coordinates and constructing the almost invariant subspace 
different from theirs, we obtain the reduced Hamiltonian which does not in¬ 
clude the vector potential terms. Using the reduced evolution we also obtain 
the asymptotic expantion of the evolution for a specific localized initial data, 
which verifies the straight motion of an atom in constatnt magnetic fields. 


1. Introduction 

In this paper we study the evolution of the system of one nucleus and some 
electrons which is neutral, i.e., the total electric charge is zero, in constant magnetic 
fields and external electric potentials. We construct the reduction scheme of the 
evolution of the system to the evolution of center of mass of the particles with 
small error terms which is regarded as Born-Oppenheimer approximation. It is 
shown that the reduced evolution of the center of mass does not depend on the 
magnetic field, which implies the straight motion of the center of mass when the 
external potential does not exists. 

The Hamiltonian of an atom with a nucleus and N electorns moving in constant 
magnetic fields has the form 

(l.I) P = —{Dxi - eiA{xi)p + - eAix^)f 

i=2 ® 

Af+1 

+ ^ ' Vij{xi — Xj) + ^ ^ Vi{xi). 
i<j i=l 

Here xi £ (resp., m) denotes the position (resp., the mass) of the nucleus, 
> 2 (resp., rUe) denote the position (resp., the mass) of electrons and Vij 
(resp., Vi) are interaction (resp., external) potentials. ei (resp., e) denotes the 
charge of the nucleus (resp., electrons) and A denotes the vector potential. 

There is a lot of literature on the quantum many body problem. Especially 
the large time asymptotics of the solution of the Schrodinger equation has been 
studied intensively and the complete classification of the asymptotic behavior called 
asyptotic completeness was proved for a large variety of potentials (see Derezihski- 
Gerard HD- We do not deal with such a problem here, but we study the phenomenon 
which happens when there are differences of the mass of the particles instead. 
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Roughly speaking the Born-Oppenheimer approximation saiys the following. Since 
the electrons are lighter than the nuclei, they move rapidly and adjust their state 
adiabatically as the nuclei moves more slowly. To develop the mathematical theory 
of the Born-Oppenheimer approximation, the ratio of electronic and nuclear mass is 
denoted by and regarded as a small parameter. Before we consider the case with 
magnetic fields, it is instructive to recall the results for the case without magnetic 
fields. In this case the Hamiltonian of some nuclei and electrons is written as 

P(h) = -h'^Ax -Ay + V{x, y), 

where we denote the coordinates of the electrons by y and that of the nuclei by x. 
Our purpose is to study the asymptotics of the solution of the Schrodinger equation 
ihdtyy = P{h)(p as h —0. From the intuitive description of the Born-Oppenheimer 
approximation above if the electrons are in bound states for the fixed nuclei, i.e. the 
bound states for Pe{x) := —Ay + V{x,y) at the initial time, we expect the electrons 
remain in the bound states even after time passes. This suggests there is an almost 
invariant subspace close to the electronic bound states under the evolution 

The almost invariant subspaces are described by the projections (see Nenciu 
mm)- If an orthogonal projection H satisfies [P(h), H] = 0{h°°), then e = 

Hg-itf _|_ 0{h°°\t\) holds which means RanH is the almost invariant subspace. 
We expect there exists such a projection H such that H — Hg = 0{h) where 
Ho = /®no {x)dx and no(x) is the spectral projection onto an arbitrarily cho¬ 
sen part of the discrete spectrum of Pe{x) separated from the other part of the 
spectrum. In Nenciu HZIIIH] the projections onto almost invariant subspaces were 
constructed by recurrence formula and another construction by pseudodifferential 
calculus was introduced by Helffer-Sjdstrand [1^ and Sjostrand m- In the case 
of Born-Oppenheimer approximation the almost invariant subspace does not seem 
to exist according to the physical intuition saying that the adiabatic decouppling 
becomes weaker and weaker when the energy increases. However for any cutoff 
function y a projection H which satisfies n]y(P) = 0{h^) is constructed by 

Sordoni [22]. Using the projection H the quantum evolution of the molecule 

is reduced to the evolution of the nuclei where G is a k x k matrix of semi- 

classical pseudodifferential operators P^(IR"), of the nuclei-variables, k 

being the rank of Hg (see Martinez-Sordoni [mils]). The symbol of G is written 
as 5 = -I- /i(x) -I- where Ik is the fc-dimentional identity matrix and 

/i(x) is a matrix of ngPe(x) in a basis of RanHg. 

Let / be a generalized coherent state of nuclei and ■)/> be a bound state of electrons. 
The initial state of the form f{x)tp{x, y) -I- 0{h) is used by Hagedorn (see Hagedorn 
[ilTlIH] and Hagedorn-Joye 0). In semiclassical limit the quantum evolution of 
generalized coherent states admits an asymptotic expantion each term of which is 
a generalized coherent state centered at the point reached by the classical flow (see 
Combescure-Robert [3] and Combescure 0)- In [16] Martinez and Sordoni applied 
this expantion to G in order to obtain the evolution like that of Hagedorn. 

We now return to the case with constant magnetic fields. Its coordinates perpen¬ 
dicular to the magnetic field stay in a bounded region as in the classical mechanics. 
When there are N particles, the center of mass of the particles does not move freely, 
so that we cannot devide the motion of the particles into the internal and external 
motion. However, if the total charge of the particles is zero, there is a subspace 
Pbound of L^(R^'^) such that internal coordinates of the particles in the state in 
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'Hbound stay in a bounded region and the particles travel to infinity across the mag¬ 
netic field (see Gerard-Laba 0). This would occur from the viewpoint of the calssi- 
cal mechanics because the Lorentz force on the particles with the opposite charge 
has the opposite direction and the particles interact with the interaction potential. 
Thus we expect when there is an atom, the electrons around the nucleus offset the 
influence of the magnetic fields. Hence it seems to be natural that the vector poten¬ 
tial terms do not appear in the reduced Hamiltonian G. The Born-Oppenheimer 
approximation with magnetic fields is also dealt with by Martinez-Sordoni |16| but 
in their construction the term h?{Dx — eA(x))^ remaines in G. This is because 
in thier coordinates where independent variables are positions of the particles, the 
effects of magnetic fields on the nuclei and the electrons are treated separately and 
therefore even if we restrict the Hamiltonian P on RanHo, the effect of magnetic 
fields does not disappear. Our purpose is to construct a reduction scheme from P 
in (HU to G without vector potential terms. 

To obtain such G we change the coordinates {xi,... ,Xn+i) to the new coordi¬ 
nates where independent variables are the center of mass x and the relative position 
y = ( 2 / 2 , ■ ,yN+i) of the electrons from the nucleus. Then the Hamiltonian P is 

transformed into a certain operator P. Next we transform P by a unitary transfor¬ 
mation V = exp (^—ieA{x) ^ VPV*. The transformed Hamiltonian 

is written in the following form: 

N+l N+1 

P = h^Dl - 4h^e A{y^)D, + 51 v)- 

i—2 i—2 

Here Li = Dy. — eA{yi), the potential V{x,y) is devided as V{x,y) = Vo{x,y) + 
V'{x,y), where Vq is the zeroth order term with respect to h, and Q is a certain 
operator on We regard Pe{x) = + K)(a;,y) as the electric 

Hamiltonian and denote the spectral projection onto the chosen part of the discrete 
spectrum of Pe{x) by no(a;). 

Our main results. Theorems 12.4112.81 in Section [2] can be illustrated as follows. 
For any cutoff function y G G^{'R.) there exists an orthogonal projection n(/i) 
onto the almost invariant subspace close to Hg, that is [P, n]x(P) = 0{h°°) and 
H — Ho = 0{h). There also exists a reducing transformation 

W : ^ (p2(M3))®fc, 

where k = RankHo, and reduced Hamiltonian G which is k x k self-adjoint matrix 
of operators on L^(]R3) without vector potential term such that the restriction U 
of W to RanH is a unitary operator satisfying 

e-^*P/>^Uipo = U*e-^^^^^UUipo + O{\t\h°°\\ipo\\), 

for (fio G Ranx(P). This means that the motions of all the particles are reduced to 
the motion of the center of mass. Moreover, if fc = 1 and the initial data is 

where H := V*nV and ui G Ran(V*no), we have the following expansion for the 
evolution 

3(A-1) 

= e*'^*/'* 5^ c^(t; 
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where St and are constants depending on t and h, is the generalized 

coherent states centered at the point reached by the classical flow of the symbol g 
of G, and Vf, € )). 

The difficulty in our construction of 11 is that when we expand the symbol of 
the resolvent (P — z)~^ in formal power series in h as z)h^ where qj 

are operators on L'^{Ry^), the power of y contained in qj{z) becomes of higher 
order as j increases. This difficulty is overcome by the exponential decay of the 
eighenfunctions of Pg and the fact that the line integral if^ := qj{z)dz where T is 
a certain curve can be written as the sum of the terms containing IIo. To prove the 
boundedness of x(P)IIj where flj := Op'^^nj), we also use the similar but slightly 
different expantion of the symbol of (P — i)~^. 

The result of Gerard-Laba [5] mentioned above also indicates that if the total 
charge is zero and the initial data is in a certain subspace of L‘^{R^), the particles 
can travel across the magnetic field. This result holds without any approximation 
but in our Born-Oppenheimer framework the localization is clearer, and although 
Gerard-Laba [5] considered the case without external potentials in which the pseu¬ 
domomentum of the center of mass commutes with the Hamiltonian, our result 
covers the case with smooth bounded external potentials. 

The content of this paper is as follows. In sec. 2 we transform the Hamiltonian 
into the form without vector potential term of nuclear variable, introduce our as¬ 
sumptions and state our main results. In sec. 3 we prove the main results. In the 
Appendix we collect the results we need on the pseudodifferential operators with 
operator valued symbols. 


2. Some preliminaries and main results 

We suppose the magnetic field is parallel to the third axis, so that the vector 
potential is written as 


/O -b 0\ 

A(x) = I ^ 0 0 j a:. 

\0 0 0 / 


where 6 > 0 is a constant. Setting the mass of electrons to 1 we introduce new 
coordinates (x, 1 / 2 , • ■ •, Un+i) = {x, y) x R^^ by setting 


x=- mxr 


yi = Xi- Xi, 


N+l \ 

-l-^Xi), M = m + N 

i=2 ) 

2 < i < IV-k 1. 


Here M is the total mass, x is the position of center of mass and yt is the position 
of the electrons relative to the nucleus. With this choice of coordinates, we have 
the Hamiltonian 


P = 


h^Di 


N+l W-l-1 

2h^e ^ ^ A(]ji)Dx + ^ ^ Li{x)'^ -I- h^Q -I- V (x, y) 

i—2 i—2 
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where Li{x) = Dy^ - eA{yi + x), and writing / = Vi^ V{x,y) and 

Q = Qi + Q 2 are as follows. 

N+l 

V(x,y) = Vii{y^) + Y 

i=2 2<z<j 

A+1 

+ Vi{x- h^f) + ^ V* (x + y* - h^f) , 

1^2 

N+l /N+l \ /N+l \ 

2 Y Lzix) Y (y^ + ) + XI ( 2 ^* + ) 

i=2 \ i=2 / \ i=2 / 


Ql — 

Q 2 = 


l-NK^ 

1 

l-Nh? 


N+l 


+ 2 X eii(a;)^(/) + h^e^A{ff. 

i^2 

We apply the unitary transformation V := exp ^—ieA(x) 2/*) 


verse to P and obtain 
p = VPV* 


N+l 


Af+1 


= h^Dl - ih'^e Y Myi)Dx + X ”*■ + '^{x,y). 

i—2 i—2 

Here, Li = Dy. — eA{yi), and Q = Qi + Q2, where 


Ql — 

Q 2 = 


1 


1 -iV/l2 
1 

1 -iV/l2 
N+l 


/N+l \ ^ 

X ’ 

^ i=2 / 

N+l /N+l 


/N+l 


^ X X (2/* + + X (2^* + 


i=2 \ i=2 


\ i=2 


+ 2 ^LiA{f) + 2h^e^A{f)^. 


i=2 

We write 

Af+l N+l 

yo{x,y) = Y yi^iVi) + X - Vj) + ^i(a^) + X + 2/*)- 

i=2 2<i<j i=2 

Note that Vo{x,y) is the zeroth order terms with respect to h of formal taylor 
expantion of V(x, y). We regard Pe{x) = Y + ^oix, y) as the electric Hamil¬ 

tonian. 


Remark 2.1. If (/? is a solution of ihdtip = P^p, then V*ip is a solution of ihdtp = P<p. 
Since V* is multiplication of the complex number with modulus 1, V*(/? has the same 
modulus as ip. 

We suppose the potentials Vij and Vi satisfy the following assumptions. 
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(HI) (i) Vij are real valued function A-bounded with relative bound smaller 
than 1. 

(ii) Vi G C°° are real valued function and for any a G there exist a 
constant Ca such that 

\d^V{r)\< 

By (HI) P is well defined as a self-adjoint operator with the domain {u G L^(]R^(^+^)) : 
{{Dxi — eiH(a;i))^ + — eA(xi))'^)u G so that P is also 

self-adjoint with some domain. Pe{x) can also be regarded as a self-adjoint operator 
with the domain {u G L^(IR^^) : ^ ^ We also suppose that 

(H2) The spectrum a{Pe{x)) is the union of two disjoint components crj(x), j = 

0,1, such that cro(x) is a part of discrete spectrum of Pe(x) with the corre¬ 
sponding subspace of L^(]R^^) being finite dimentional and there exists a 
number d > 0 such that 

inf dist(CTo(a::), tTi(a:)) > d. 

Remark 2.2. By (H2) it is easy to see that there exist a continuous family of loops 
7 (a;) which encloses croix) for each x and inf^-gRa dist( 7 (a;), (T(Pe(a;))) > d' for some 
d' > 0. 


We denote by no(x) the spectral projection of Pe{x) corresponding to CTq. We 
also suppose the following assumption. 

(H3) Ranno(x) is spaned by a orthonormal basis {ui{x, y),..., Uk(x, y)) G Ly) 

such that 


/R3N 


|u,(x,y)|2e2“IJ'ldy <C, 


where constants C > 0 and a > 0 do not depend on x. 


Remark 2.3. If Vi{r) =0, 1<*<7V-|-1, then ai{x) and no(a::) don’t depend on 
X, so that the assumptions (H2) and (H3) are obviously satisfied. If the fluctuation 
of Vi is small enough, by the upper semicontinuity of the spectrum the eigenvalues 
in (Jo do not cross each other, so that the corresponding eigenfunctions are smooth 
orthonormal basis. Moreover if Vi are periodic, the eigenfunctions decay uniformly 
with respect to x since fie in the proof of Agmon [1] Theorem 4.1] can be chosen 
uniformly in x. 


Our main results are concerned with the almost invariant subspace which is close 
to electronic eigenspace. In the following theorem, a = 0{h°°) means a = 0{h^) 
for any A G N. 

Theorem 2.4. Assume (H1)-(H3) hold true. Then for any $ G (^“(R) there 
exists a orthogonal projection n(h) on L^(R^(^+^)) such that 

l|n - no||_f^2(R3(N+i)) = o{h), 

and, for any x € Co°(I^) such that x4> = Xj have 

||x(p)[n,p]|L2(R3(.+o) + ll[n,P]x(P)IL2(R3(.+o) = o{h^). 

Theorem 2.5. If ip is the solution of ihdtp) = Pip with initial data ipo satisfying 
xiP)ipo = To /o*' some x G (^“(R) such that x<P = X> then 

(2.1) ip = e-^^P^/^Uipo + - n) 7 .o + O(|t|h“||^o||), 
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where Pi := IlPn and P 2 = (1 — n)P(l — 11) is self-adjoint on a domain containing 
D{P). 

We can reduce the evolution to an evolution on the space of only 

nuclear variables. To state the next result we use h-admissible operators. For the 
definition of h-admissible operators see the the Appendix. 

Theorem 2.6. There exists a h-admissible operator 

W : ^ (p2(R3))efe 

with operator valued symbol and kxk self-adjoint matrix G of h-admissible operators 
on such that the restriction U ofW to Ranll; 

U : Rann ^ 

is a unitary operator which satisfies 

uPiU = Gun, 

so that = U*e~^*^^^UU. The symbol g{x, of G has the following form: 

gix,f) = + fxix) + ^h^ gjix,f), 

i>l 

where pL(x) is the matrix 0 /no(x)Pe(a;) in {ui{x),... ,Uk{x)) and g{x,ff) does not 
contain vector potential terms. 

Remark 2.7. When more than one nucleus exists, these results with the operator 
G without vector potential term does not seem to hold because each electron is 
possessed by more than one nucleus at the same time and the electrons around 
each nucleus does not necessarily offset the influence of magnetic fields. 


The next result is concerned with more explicit expression of the solution of 
the Schrodinger equation for a special initial data. Let us write Uj{x) := V*Uj{x). 
Then ({(i(x ),... ,Uk(x)) are the orthnormal basis of Ranno(x) where 


no(a:) := V*no(a:)V. 

which is the spectral projection of Pe{x) = Li{xY + Voix,y) corresponding 

to ao(x). 

Let at = {xt,f,t) be the solution of 


( 2 . 2 ) 


xt = ^{xt,^t), ft = 


starting from initial data oq = (xojfo)- Let (jin, Cn), n = 1, 2,3 be the independent 
solutions of 


(2.3) 

with initial data 


P =JMt 


{rin)j\t—0 — {fn')j\t—0 — ihjn^ 


where {rin)j is the jth component of 


J = 


0 I 
-I 0/ ’ 
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I being the unit matrix and Mt is the Hessian of g at at'. 


(Mtkj 



oc—oct 


By Theorem 12.51 and Theorem 12.61 as in Martinez-Sordoni [161 Theorem 11.3] (see 
also Combescure-Robert m), we have the following theorem. 


Theorem 2.8. Let k = 1 and (po £ be as follows 

where H = V*nV, and x = 1 near + g{xo). Then there exists C > 0 such that 
for any integer J > 1 one has 

3(J-1) 

(2.4) X! Cti {^\ h )( j )^^ tv ^,{ x ) + Olh -’/^), 

/i—0 


where 


Jfj. 

i=o 

c^j are polynomials with respect to d'^g{xt, ^t) and Kepn, Inir?„, Re(C„, ImC„ 1 < n < 

3,dt := J*{xs^s - g{xs,^s))ds + {xofo - Xt^t)/2 andv^ £ (^“(M^; L^(M^^)). The 
estimate is uniform with respect to {t, h) such that h > 0 is small enough and t < 

C~^\nj^. Moreover, the lowest order term with respect to h is ^‘^^ui(x)). 


3. Proof of the main results 

The outline of the construction of H can be illustrated as follows. First we obtain 
the expression of the symbol of the resolvent {P — z)~^ of P as formal power siries 
tlj{x, C) z)h^ with respect to h. Secondly we integrate each qj with respect to z 
along the loop r(a;, enclosing the set {z : z—^ £ cro(a;)} and denote it by k^{x, ^). 
Thirdly we quantize kj and multiply the cutoff function 4>(P) to make each term 
bounded. Fourthly we make a resummation and symmetrize the operator to obtain 
a operator H^. Finally integrating the resolvent of H^ along the loop enclosing the 
point z = 1, we obtain H. 

We denote by p{x,f', h) = ^ MVi)^ + h'^Q + ^ (x, y) 

the symbol of P and by po{x, := + ES^ + ^o{x, y) its principal symbol. 

Let 'y{x) be the continuous loop as in Remark Let us consider H := {{x, z) £ 
K® xC ; z— S 7 ( 2 ^)}- For (x, z) G Q, po{x,f) — z is invertible and qo := {po — z)~^ 
is smooth and bounded. We define the symbol 

r{x,^',h,z) = ^P{x,f',z)h\ 
i>i 

as in Sordoni [^ (see also Nenciu-Sordoni [Il)by 

(p(x, h) - z)#qo{x,^', z) = 1 - r(x, h, z), 


aix,^)#b{x,f) = 


h\o‘+h\(—\\\o‘\ „ 


where 
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at a formal series level. We define 


q{x,^;h,z) = ^q>{x,£,\z)h^, 
i>i 


by 


q := qQ + qQ#^r*^. 

j>i 


Here where the number of ^ is j. One can check that, for j > 1, 

qj is given by the sum of the terms of the following form 


(3.1) 

i=l 


where C is a constant, 1 < m < 2j, = \Y^T=i Pi\ — j Si is one of the 

following 


(3.2) 


V+l /N+l \ 2 /N+1 \ ■ 

C, |a| = 1 , L,A{yi), ( X! ’ [Yl ^MVi) ) 

i=2 \ i=2 / \ i=2 / 

( A+1 \ ^ N+l /N+l \ ^ 

Yy^j +Y +j/*) I ^2/^I , 7€N^. 

i=2 / i=2 \ i=2 ) 


Since the operators do not commutate with each other, the order of the product 
is specified as n+Li((^“'^f'sO<?o) = (9“i5^^si)(7o • • • (9“™9^’"Sm)<?o- Let us set 

r(a;, 0 := {2: e C; z - e 7(2;)}, 


and let us define 

(3.3) nj{x,^) = -^(p qj{x,^;z)dz. 

27!" Jrix,i) 

Then by the change of the variable z ^ z — we see 7ro(x, ^) = no(x). 

If a potential V in K.^ is A-bounded with relative bound e, then V{yi) and 

V{yi — yj) are relatively bounded with respect to X)i =2 Lf with relative bound e. 

Thus we see easily the following lemma. 

Lemma 3.1. Let Re{x; z) = {Pe{x) — z)~^ he the resolvent of Pe- Then 

(i) {Li)kRe{x;z) and LiLjRe{x] z) are hounded as operators on uni¬ 

formly with respect to x where {Li)^, k = 1, 2, 3 denotes the kth component 
of L,. 

(ii) {Li)kqo, LiLjqo, ^^qo and \^\{Li)kqo are hounded as operators on 
uniformly with respect to x,^. 


To quantize itj we need to prove that they satisfy the condition of the symbols 
(see the Appendix). 

Lemma 3.2. For any j GN and r, s € R, kj{x,f) G S'2(IR^ x £(Ly’’', Ly’’'+®)) 
where := {u : (1 + \y\^Y^'^u G L^}. 



10 


SOHEI ASHIDA 


Proof. By (13.ip qj can be written as the sum of the terms of the following form 

m 

(3.4) 

where C is a constant, G /3 G and aiRe{x;i) G uniformly 

with respect to x. It is easy to see \f3\ < j from the construction of qj. By the 
change of the variable z ^ z — ^qo YYiLi{y°‘'<^iQo)dz becomes 


P m 

(j Re{x\z)W{y°‘'aiRe{x]z))dz. 


’l(x) 

Let 7 i(a;) be a slightly larger loop enclosing ^{x) and Zi is a point on 71 (x). By 
the resolvent equation we have 


® ® Re{x;zi)Re{x;z) aiRe{x; z))dzdzi 

J^{x) Jjl(x) 

= ® ® — - — {Re{x-,z)-Re{x-,zi))T\{y°‘'aiRe{x;z))dzdzi. 

j~f(x) j^-[(x) ^ ^1 


h(x) Jll(x) 

Since 71 (x) lies outside 7 (x), we have 


f f —— Re{x;z)Y\{y°‘'aiRe{x;z))dzdzi 


« m 

= -2-ni ® i?e(x; z) J|(y“*aii?e(a:; z))dz. 


Thus we have 


® ® Re{x; zi)Re{x; z) aiRe{x] z))dzdzi 

d-t{x} J-fi{x} .■_n 


P m 

= -2Tri ® i?e(x; z) J^(?/“’aii?e(x; z))dz. 

Jj{x) 


i=l 


- ® ® —-—i?e(x;2i)TT(i/“‘aii?e(x;2;))dzdzi. 

J l{x) J '-ii(x) Z' — Zl 

Transposing the terms and using no(x) = ^ — z)~^dz, we have 

P m 

® Re{x;z)Y\_iy°‘'aiRe{x;z))dz 
dxix) 

P m 

= ® no(x)i?e(x;z)]^(y“*aii?e(x;z))dz 

Ji{x) 

i 


- Re{x; zi) ]T(?/“‘aii?e(x; z))dzdzi. 

■yix) djiix) Z — Zi 


i=l 
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Repeating the same procedure for the rest of Re{x; z) we have 

P m 

® Re{x;z)Y\_(y°‘*aiRe{x;z))dz 

\2 t / Jj(x)J'yi{x) / 


(3.5) 


■no(a;) Re{x-,z){y°‘'ai)\Re{x]z)dzdzi---dzi 

\i=f+i / 

m+1 


m+1 

n 


1 


7(2:) ./7i(a;) •t7m+l(a:) “1 ^ 

m 

n<« e{x] Zi){y'^'ai))Re{x] Zm+l)dzdzi ■ ■ ■ dZjn+1, 


where ji{x) are the same loop as 71 (a;). Since 7 i(a;) lies outside 7 ( 2 ;) the integration 
with respect to z in the last term of the right hand side vanishes. 

By the assumption (H3) y^Ilo{x) is a bounded operator on uniformly 

for X. By 


[{yi)k,Re{x; z)] = Re{x; z)[Pe, {yi)k]Re{x-, z) = -2iRe{x; z){Li)kRe{x; z), 
the remaining term of (|3.5I1 are the sum of the terms of the following form; 


(3.6) 




7i(2 


Rj^ln 




y^Uo{x)y^ Re{x; z)ci J Re{x; z)dzdzi ■ ■ ■ dzi, 


where C is a constant, y^Ilo{x)y^ is bounded uniformly in x and Re{x]i)bnj, 
Re{x]i)cj are bounded uniformly in x. Since 


{Pe{x) - i)Re{x; z) = 1 + (z - i)Re{x; z), 

is uniformly bounded for x, z € 7 (x) and (z — Zk)~^ is uniformly bounded in 
z € 7 ( 2 :), Zfc G Jkix), (13.61) is bounded uniformly in x. Thus by (13.4|) || 71 ^( 2 :, C)IUg < 
(7(1 + 1^1)^ uniformly for x, where C > 0 is a positive constant. In the same way 
we can see for all a £ and 13, S £ there exists a constant Cc,,y,s such 
that \\y°'d!^d^TTj{x,^)\\L 2 < • This completes the proof of the 

lemma. □ 


Remark 3.3. Since in Martinez-Sordoni m they use the coordinates where inde¬ 
pendent variables are positions of the particles, the power of y appears only in the 
electric Hamiltonian and is not included in qj. On the other hand in our coordinate 
the power of y is included in A{yi)Dx, Q, and the taylor expantion of 

Vi- Thus in our case the power of y in qj becomes of higher order as j increases, so 
that the boundedness of ifj is not obvious. 

We define 

lij = OpI{tTj{x,0)- 
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The operator 11^ is not bounded but by a localization in energy we obtain a bounded 
operator. 

Lemma 3.4. For any 4) G C^(K.), <l>(P)nj is bounded in L'^. 

Proof. We can write 

(3.7) $(P)n, = <i>(P)(p - iy{p - 

where 4'(P)(P — iy is bounded. We set 

p\x,^) :=po(x,0 + h^Qi- 

We define 

%{x,£,]h) := (,p'{x,C) - 

and r'{x, h) by 

{p{x, f-,h)- h) = l- r'{x, h). 

Since 

= -qQ{x,f,-,h)d^^VoqQ{x,^-h), 

and 

= -Qo{x,^;h)^kqo{x,£,\h), 

by Lemma [SHI we see r'{x,f;h) G x £(Ly’®, for any s G M. 

Since 


/j-i ^ 

{l-r')#ij2r'*^ 


= I — r‘ 






we have 


{P - ^)Opy 0 # r'*^ ^=l-Opl {r'*i) . 

Multiplying the both sides of the equation by (P — i)~^ from left we obtain 

(P _ ,)-i = (P _ ,)-iop- (P#.) + Opy (^q'oix,0# ) ■ 

Hence we have 
(3.8) 

{P-i)-%=Y.{P-^)-^+'^Opl {r'*^) iopl U{x,i)if (E^'n ) 
fc=0 V V \i=0 /// 


fi,+ (opy (^q'o{x,0# 


n,. 


Since for any r G 


and 
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by Lemma l3.2l we have 




0-1 






\i=0 


Thus by Theorem lA.41 it follows that 


^i-i 


Opl {r'*^) Op't: q',{x,0# Y^'*^ 


0=0 


n, G£(L^(R^;L^(Ry))) 


^i-i 


(3.9) 


Op)r 9 o(a;,e)# [Y^'*' n, € £(L 2 (M 3 (iv+i))), 


0=0 


By (1^ . (15^ and (1531) . we have G £(L^(]R^(^+^))) 


□ 


Remark 3.5. In Martinez-Sordoni m they considered only localized initial data, 
so that unboundedness of the vector potential term is not important. On the other 
hand, in our case we do not restrict the initial data and therefore, because of the 
existence of the vector potential terms in P in cu, Dx is not P-bounded. Thus, 
to prove (P — is bounded, we need to expand the resolvent at the expense 

of the presence of the power of y. Since tTj includes no(a;), the power of y does not 
prevent (P — from being bounded. 


To construct 11 we make a resummation. Let p G C))°(]R) be a function such that 
0 < p < 1, supp p C [—2, 2] and p = 1 on [—1,1]. It is easy to see that Pnj$(P) 
and ^{P)fljP are bounded in P^(]R3 (^+i)) a^d there exists a decreasing sequence 
of positive numbers (ej)jgN converging to zero such that for any j G N 

(i -p (D) (||n,$(p)|| + ||$(p)fi,|l + ||pn,4>(P)|| + ||<i>(P)n,p||) < h-\ 

Let us define 

n<i>(p) := no$(P) + E (i - (I)) 

t>i 


$(p)n := $(P)no + Y{^-p{^)) 

3>i ^ 

Then it is easy to prove that for any iV G N there exists Cm > 0 such that 


(3.10) 
and 

(3.11) 


N 


n$(p)-^hm,$(P) 

i-0 


N 


4>(p)ri-E^'‘&(^)nj 

3=0 


< CNh^+\ 


N 


pn4>(p) - E h^pfijHP) 

3=0 


N 


4 >(p)np - E h^HP)^jP 

3=0 
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where we understand IIo = Hq. We define 

n$ := $(P)fl + (1 - $(P))fl$(P) + (1 - $(P))no(l - $(P)). 

Since formally the symbol of {P — z)~^, by (13.101) and (13.111) we have 

the following lemma as in Sordoni |22) . 

Lemma 3.6. Let x S C'^(R) satisfy x4‘ = X- Then for all N gN 

||x(P)[n$,P]|| + \\[u^,P]xiP)\\ = oih^), 

||x(p)(ri| - n$)|| + ||(n| - n3,)x(P)ll = oih^). 

Since n$ —Hq = 0{h), for sufficiently small h the spectrum of 11$ is concentrated 
near 0 and 1, so that the set {z S C; jz — 1| = 1/2} is in the resolvent set of 11$ for 
sufficiently small h. We define 

11:=® (n$ — z)~^dz. 

J\z-l\ = l/2 

Then as in Sordoni [22] we have 

n-n$ = ;:^(n|-n$) / (n$-z)-i( 2 n$-i)(i-n$-z)-i(i-z)-idz. 

27r J\z-i\=i/2 

Theorem 12.41 follows from this formula and Lemma 13.61 


Proof of Theorem \2.5\ First, we prove that IIPII and (1 — n)P(l — 11) are self- 
adjoint on a domain containing D{P). Since it can be proved that Pn(P — i)~^ 
is bounded as in the proof of Lemma [3.41 we have nP(P) G D{P). Hence HPH 
and (1 — n)P(l — n) are defined on D{P). Since P is lower semibounded, P is 
also lower semibounded. Thus HPH and (1 — n)P(l — H) are lower semibounded 
so that there exist Friedrichs extentions of HPH and (1 — n)P(l — H) with their 
domains containing D{P). 

Multiplying both side of ihdt<f = Ptf by H we have 

ihdtLly} = PiHt/? -I- Piv?, 

where Pi = Hp, P]x(P) = 0{h°°). This equation can be re-written as, 

and thus integrating from 0 to t we obtain, 

+ O(|i|/i°°|po||)- 


In the same way we obtain 

(1 - H):^ = - n)(^o + O{\t\h^\\^o\\). 

Combining these two equations we obtain dm). 


□ 


Proof of Theorem \2.6\ Since H — Hq = 0{h), for h small enough the operator 
U := (HoH + (1 - no)(l - n)(l - (Ho - 

can be defined. W is a unitary operator and it maps RanH onto RanHn (see Kato 
[HI Chap.1.4). 
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We can write $(P)nj = ^{P){P — iy{P — By the proof of Lemma [3^ 

{P — is a h-admissible operator with a symbol in x £(Ly, 

for any s € K. Writing $'(a:) := $(a;)(a: — iy we have 

(3.12) <i>'(P) = ^ ^ d^jdizyp - z)-^dzdz, 

where $'(z) is an almost analytic expantion (see e.g., Martinez El)- Expanding 
{P — z)~^ in (|3.12p as in the proof of Lemma 13.41 we can see that $'(P) is h- 
admissible, so that n$ is a h-admissible operator on By the equation 

n = ® {z — n$)(iz 

d|2-l| = l/2 

p oo 

= f {z- no)-i V((n3. - Uo){z - Uo)-y’^dz, 

J\z-l\ = l/2 

we see that U is a. h-admissible operator on L'?,) and its princial symbol is 1. 

We define W by 

wy = {Uy,ui{x,y))L2(^3N) © ... 0 [Uy,Uk{x,y))L2(s3N). 

where (ui,..., Uk) is the orthnormal basis of Ranllo in (H3). Since W*(ai © • • • © 
ak) = U*{aiui + • • • + akUk) for Ui G L^(]R^) we obtain 

w*w = whqU = n, ww* = 1 , 

which implies the unitarity of the restriction {7 of W to Ranll. 

Difining G := UPiU* it is easy to see that G is a h-admissible operator on 
Ly) and its symbol have the following form: 

5(a;,0 = ^^4 +/i(a;) + ^h^gj(a;,0- 

i>l 

We shall prove G is self-adjoint with the domain ?7(Rann fl D{Pi)). For this 
purpose we only need to prove Pi is self-adjoint with the domain Ranll n P(Pi) as 
an operator on Ranll since U is unitary on Ranll. Let u,w G Ranll and suppose 
{u,Piv) = {w,v) for all v G Ranll n D{Pi). Then for any v G D{Pi), (u, Pi flu) = 
(u’,!!'!;). Moreover since Pi = flPII and w G Ranll, we have (u, Pi(l — !!){;) = 0 
and {w, (1 — n){;) = 0. Thus we have {u,Piv) = {w,v) for any v G D{Pi). Since 
Pi is self-adjoint, we see u G P(Pi) and therefore u G Ranll fl D{Pi). Thus Pi is 
self-adjoint with the domain Ranll n P(Pi) as an operator on Ranll. 

Since G is self-adjoint we can define and obtain 

This establishes the theorem. □ 

Proof of Theorem \2.8i In the same way as Martinez-Sordoni [THl Theorem 11.3], 
for ifo = (7rh)-44nx(P)(e“f“4-(^-^o)V2/iy^(^)) we obtain the following: 

3(J-1) 

(3.13) e-*‘^/Vo = ^ c^(t; h)</.^,*u^(:r) + 

fi—O 

where 5t,c^{t\h) and 4>^^t are satisfying the conditions in Theorem 12.81 and v^{x) 
is satisfying the same condition as v^{x). 
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Since V*</Jo = <^o, multiplying both sides of (|3.13p by V* we have (I2.4p with 
Vfj,(x) = V*v^(x). Since Vfj,{x) contains ui, we have 

d,^{V*v^{x)) = -ie (a yn V*v^{x)+V*d,^v^{x) e C°iRl;Ll). 

\ i=2 / j 

In the same way V*i;^(x) is differentiable any number of times which completes the 
proof. □ 

Appendix A. Pseudodifferential operators with operator valued 

SYMBOLS 

We introduce the classes of operator valued symbols we use. 

Definition A.l. Let "Hi and 'H 2 be Hilbert spaces and Hq a positive constant. A 
function a{x,^;h) G x x {0,ho]-C{'Hi,'H 2 )) is said to be in x 

C{'Hi,'H 2 )) if for any a G N" and /3 G N™ one has 

(A.l) sup < 00 . 

(x,4,/t)GR"-+^x(0,/to] 

Definition A.2. For a{x,y,^) G x K"; "^ 2 )) and u G S{RA\'Hi) we 

set 

OPh{a)u{x;h) = J d£,. 

Op™ (a) is called the semicllasical pseudodifferential operator with operator valued 
symbol. 

We can see that Op™(a) is continuous 5(R";Hi) —5> S(R™;'H 2 ) as in Martinez 
[l4] . Since the formal adojoint of Op™{a) is Op™{a*), we can extend it uniquely to 
a linear continuous operator S'{R™]'Hi) —J- 5'(R";7t2). 

For the composition of pseudodifferential operators we have the following theo¬ 
rem as in the scalar case (see Hormander [11)1. 

Theorem A.3. //oi G S'™! (M” x R”; H 2 )), 02 G S™=!(K" x R"; £(^ 2 , ^ 3 )), 

then as operators onS(R";7ti) orS'(M";'Hi) 

Op)('(a2)OK(ai) = 

where b G 5 '™i+"! 2 (K" x R”; £( 7 ^ 1 , Tta)) is given by 

6(x,C) = =: ai#a2, 

v=i 

and has the expansion 

u\ci+P\( _idal 

(o-\\a+0\ \R\ {9'^d^a-2{x,0){d1ddax{x,i))+h^+^rN{x,C). 

|a+/3|<7V 

where rN G ^ R". 

The L^-boundedness can be established as in the scalar case. 

Theorem A.4. Let a G S°(R"' x R"';£(7t)). Then Op™{a) is continuous on 
L^im™;n) and 

l|0£“(a)llc(L2(R";W)) <Cn{ ^ ll^“^fa||L“(R"xR";£('H)) 

\ |a+^|<M„ 
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where the positive constants Cn and Mn depend only on n. 

We also use the useful notion of h-admissible operators (See Martinez-Sordoni 
[l6] Appendix). 

Definition A.5. Let m € M and let Hi and H 2 be Hilbert spaces. An operator 
A = A{h) : —>• with h € (0, ho] is called h-admissible (of 

degree m) if, for any A > 1 

N 

= X! th)) + h^Ajv(h), 

j=o 

where Rn is uniformly bounded from 7J'"(R";?^i) to L^(M";H 2 ) for h € (0, ho], 
and, for all h > 0 small enough, Oj G S'"‘(M^"; £(?ti, "^ 2 )). 
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